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Abstract 

We consider a Gaussian random matrix theory in the presence of an exter- 
nal matrix source. This matrix model, after duality (a simple version of the 
closed/open string duality), yields a generalized Kontsevich model through 
an appropriate tuning of the external source. The n-point correlation func- 
tions of this theory are shown to provide the intersection numbers of the 
moduli space of curves with a p-spin structure, n marked points and top 
Chern class. This sheds some light on Witten's conjecture on the relation- 
ship with the pth-KdV equation. 
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1 Introduction 



The intersection theory for the moduli space of curves with marked points has 
been studied by various methods with the aim of proving Witten's conjectures 
[1] on the relationship between the structure of this space with KdV flows. 
In particular this has led to a solution through Kontsevich's Airy matrix 
model[2]. The generalizations of this theory have been also considered also 
by Witten, who conjectured that the intersection numbers of the moduli 
space with top Chern class is described by the pth-KdV equation [3]. 

This generalization is related to the twisted N=2 supersymmetric confor- 
mal field theory, and particulary it has a close relation to the level k = p — 2 
WZW model [4]. 

In previous articles [5, 6], we have found that a random matrix theory with 
an external source gives an alternative method for obtaining the intersection 
theory of the moduli space of Riemann surfaces with a p-spin structure. 

We have discussed it through an N-k duality in the expectation values 
of the product of k characteristic polynomials of N x N matrices . After 
tuning appropriately the external source to a critical value, we obtain the 
intersection numbers of the moduli space from the Fourier transforms of the 
correlation functions. We had studied in the past random matrix theory with 
an external source as a way to tune higher edge singularities, including gap 
closing cases [7, 8]. Through this duality, plus tuning of the source, we show 
here that one obtains the generalized Kontsevich models. 

In a previous article, we have computed the intersection numbers with 
one marked point for the generalized Kontsevich model with p-spin curves by 
the replica method applied to a Gaussian random matrix theory. We have 
also found, by an alternative method, that the Fourier transform of the one 
point correlation function U(s) provides those intersection numbers and we 
have found an explicit agreement of the two approaches in the case of one- 
marked point and arbitrary genus when p = 3 [6]. This was checking that 
a Gaussian random matrix model with an external source is a dual (mirror) 
model to the generalized Kontsevich models of p-spin structures with one 
marked point. (We will also show below that it is true for p > 3.) 

For p=2, which is Kontsevich original Airy matrix model, the intersection 
numbers have been computed explicitly, for a few marked points, from a 
Fourier transform of the correlation functions at the edge of the spectrum of 
the density of eigenvalues, in a random matrix model by Okounkov [9, 10] 
and in our previous article [5]. Here we extend the analysis to the case of 
p ^ 2, again from the edge singularities of the correlation functions in a 
matrix theory. 

The intersection numbers are characterized by the genus g and marked 
points n with a "spin-structure" labelled by a non- negative integer j, (j = 



1 



0, 1, ...,p — 1). Higher Airy matrix models for curves with a p-spin structure 
are given by [2] 



Z = — [ dBexp[—^—tr(B p+1 - A p+1 ) - tr(B - A)A P ] (1) 
Z J p + 1 

where 

r p ~ x 1 

Z = dBexp^tr- A j BA p -^ l B] (2) 

J j=0 2 

The free energy, the logarithm of the partition function Z, is the gener- 
ating function for the intersection numbers < l\r m j >, 

^=£<n^->im (3) 

d m ,j m,j m,j m >3- 

where 

j-p-m(p + 2) n l— w ^ 1 

tmj = {-p) ^ n ('P + 3 + !) tr x^w ( 4 ) 

According to Witten [3] the intersection numbers is given by the top Chern 
class Ct and by the first Chern class C\ as 

If n 
< r mujl ■ ■■T mnJn >= — I c T (ji, ■ ■ -,j n ) II C i(^) mfc ( 5 ) 

with the condition which relates, for given p, the indices to the genus g of 
the surface, 

n 

(p + 1)(20- 2 + n) = £(prni +£ + !)■ (6) 

i=i 

For more precise definitions and recent studies of the intersection num- 
bers, we refer to the literatures [17, 18, 19, 20, 21]. 

The higher Airy matrix model of (1) is known to corresponds to a {p — 
1) matrix model, which is conjectured to satisfy pth KdV equation. The 
Virasoro equations for this case are given in [22]. They have been used for 
computing recursively these intersection numbers ; here these numbers are 
given directly through a generating function. 



2 One point correlation function 

In a previous article [5] we have used explicit integral representations for the 
correlation functions [12, 13, 14] for a Gaussian unitary ensemble (GUE) of 
random matrices M in the presence of an external matrix source. 
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The probabililty distribution for N x N Hermitian matrices is 

P A ( M ) = -L e"T ^-NtrMA (7) 

Za 

We denote the eigenvalues of the external source matrix by a a . 
The Fourier transform of the one point correlation function, the density 
of state p(X), is 

U(s) = / — e lXs < tr£(A - M) > (8) 
J 2ix 

where the average < ■ • • > is taken with respect to the probability distri- 
bution Pa(M). An exact integral representation is known [12, 13, 14] for 
arbitrary a,j : 

U(s) = — ef s2 I — ft (1 - ^—)e Nsu (9) 

In the following we shall consider special cases in which the degeneracies 
of the distinct eigenvalues are proportional to N. In particular we consider 
(p — 1) different eigenvalues a a , (a = 1, ...,p — 1), each a a being times 
degenerate. Then above expression reduces then to 



We now consider the large-N limit, in a regime in which s and u are of 
order N~( p+1 \ The distinct eigenvalues a a are taken to be all of order one. 
Let us expand the integrand 

1 n s 2 r du N ^ ~ Au + s) m u m . Ar . 
U(s) = —ei s & exp > > (- ) + Nsu] (11) 

In the large-N limit the summation over m may be truncated at m < p + 1 
since higher orders vanish in that limit. We now specify the (p — 1) distinct 
eigenvalues of the external source by the {p — 1) conditions : 

£i=P-l, E^ = °> (m = 3,4,...,p) 
£ 4t ^ 0. (12) 
With these conditions, in the regime of interest, one has 

= £)] / ^expC K « - (« + (13) 
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with 

N p ~ l 1 

The integration over u takes a more symmetric form after the shift u — > m — |s 

[/(«) = J- e -^E^i £ / *i e -o/(«) (15) 

v ; iVs 7 2vri K J 

with 

/(«) = (« +1^-^-1^ 

= E( p ^ 1 )£(i-(-i) m )^ 1 - m ( 16 ) 

where the interval of the integration of u is (— ioo, +ioo). 



3 Edge singularities 

In the large N limit, the density of state p(A) has a finite support, and thus 
develops singularities at the edge of the distribution. We consider now the 
nature of this singularity for the p-spin case. 

In the large N-limit the Green function G(z) = ^(tr^j) is given by a 
simple equation due to Pastur [15] (see also [16]) 

For a source which consists of (p — 1) distinct eigenvalues, each of them 
degenerate N/ (p — 1) times, this reads 



I 

'G(z) 



GW =^iS^^ (18) 



If we expand this resolvent in powers of — , and use the conditions (12), we 
find 

i p-i i i p-i i 

" A(E^h)(--^ (19) 
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in which we have neglected terms of order (z — G) p+1 or higher. From the 
first condition (12), we have 

^-A(Ei)(^'-GV)) (20) 

P 1 a =l a a 

where z' = z + -Xr £ — and G'(z') = G(z) + -srr- The Green function 

p-1 ^ a a V / V / p-1 ^ a P +1 

has a singularity proportional to z n / p ; thus the density of state p(X) = 
— ^lmG(A) has an edge singularity characterized by an exponent K 

Note that the conditions (12) for the a a admit several solutions. However, 
all these different choices lead to the same singular behavior. In other words, 
the singularity exponent is independent of the location of the critical points 

For instance, in the case p — 3, the (complex)solutions are 

(ai,a 2 ) = ±(1, -1) 

and 

, . , , V3 1 . V3 1 

(oi ' tt2) = ±( V! + V5 ! 'V5-V5 !) 

The dual higher Airy matrix model in (1) is obtained from expectation 
values of characteristic polynomials [5, 6], leading to 

p— 1 p— 1 • j3 P— 1 

< J] det(a a -tB)^ >=< exp[^trlog(l-— )+iV^log(n a a )\ > (21) 

a=l a=l a " a=l 

Expanding the logarithm, with the conditions (12), we obtain in the large N 
limit, the higher Airy matrix model (1) as explained in [5, 6] (up to a change 
of the normalization factor). 

4 Intersection numbers of one marked point 
for £>-spin curves 

We now derive the intersection numbers of one marked point from the asymp- 
totic series expansion of U(s). 

(i) p=2 We begin with the simple edge of the semi-circle law. In the 
large N-limit, in the range in which s is of order N^ 1 ^ 3 , we have 



U(s) = —e 4 s / —e 3Csu 
v 1 Ns J-oo 2m 



Ns V SCs 



K < (22) 
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where C = — y. By the change of the normalization due to the higher Airy 
matrix model of (1), we have s 3 /24 = —Ns 3 /12, and we get 



1 / 7T S 9 

1 T= (23) 

Noting that s is a conjugate variable to A, this yields a series expansion 
in inverse powers of A. From the definition (3), we obtain the intersection 
numbers for one marked point, 

< T3 »« >•= pTw' <24) 

which agrees with our previous result [6] based on the replica method. 

(ii) p=3 The critical point now corresponds to a density of states whose 
support consists of two disconnected segments, in the limit in which the gap 
closes. The intersection numbers with one marked point for p=3 have been 
obtained, for arbitrary genus g, in our previous replica article [6] : 

i rf^i 1 ) 

< T "^> ' (iwF(W (25) 

where the spin-index is j = for g = 3m + 1 and j — 1 for g — 3m 
(m=l,2,3,...). 

Near the edge point, z c = — ^ J2 — , we find in the scaling region from 
(15) 

= w 4 «) < 26 > 

where ( = —N 2 ^ 3 (A ■ 3 1 ^ 3 )~ 1 s 8 / 3 which was already established in [6]. From 
the standard asymptotic expansion of Airy functions Ai(Q, we have two 
distinct series, 

+^'(0)(C + §C 4 + 2 -^C 7 + 2 -^C 10 + •••)] (27) 

where Ai(0) = 3- 2 / 3 /r(|) and Ai'(0) = -3- 1 / 3 /r(|). 

The first series in (27) gives the intersection numbers for j = 1. From 
this series , noting that s ~ we have 

11 14 

< T6,l > 9 =3- • -, < 714,1 >g=6~ ' -, 

1 1-4-7 

< r 22 ,i > g =9= ( 12 y Q \ ' 33 ' ( 28 ) 
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From the second series, we have 

1 1 2 

< Ti,o >g=i- < T9 '° > f= 4_ (12)441 ' 3' 

12-5 1 2-5-8 

< Ti7,0 > 9 =7- ^ 77! ' -p-, < 7-25,0 >g=10- ( 12 )10 10 ! ' 3 3 -( 29 ) 

These values agree with (25) and with the results obtained by completely 
different mothods [19]. 
(iii) p = 4 

In this case, the critical values of the a a , which satisfy the conditions 
(12), are obtained as the zeros of a cubic equation, 

b a = — (a = 1,2,3) 
d 



a 

3 1 ajl 



with 



(x - b x ){x - b 2 )(x - 63) = x A + (3x 2 + + 5 = (30) 

P 2 = 9±3V6, 7 = ^ 2 -3), 5= (31) 

There is one real solution, and two complex conjugate solutions for the a a . 
(Although an analytic expression for the a a , solutions of the cubic equation 
(30) exists , we give here the numerical values). There are two classes: 

(a 1? a 2 ,a 3 ) = (±0.52523, ±0.41127 ± 0.46403i, ±0.41127 =f 0.46403i), 
(ai, a 2 , a 3 ) = (±1.0076, ^0. 71801 ± 0.33908i, ^0.71801 q= 0.33908i) (32) 

In both cases, the density of state has one critical edge, at which it behaves 
as p(A) ~ A 1 / 4 . Note that , contrary to the critical gap closing p=3 case, for 
p=4, there is just one single edge, similar to the p=2 case. 
We have 

U(s) = — e s*(2.. a ) jri5 iC / exp [-C(-v A - -v 2 )] 
iVs J -00 2tt 3 b 

Ns Jo 2tt v ; 

where C = iVX)a=i -V- In a series expansion in powers of s, we obtain 



K 4Ns JK 3 J 

X[r( 4 )(1 + 2!-4 ( l2- )+ 4T ( T2- ) V + '-- ) 

B ,iVC\i„,3 w 13s 5 , NC, s w 4:-7NC, 2 

+ ^-(^) + ^(^) 2 + 



(34) 
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For this p = 4 case, there is the overall factor e~4 2 -is s . It must also be 
expanded and combined with the series to find the relevant terms in the s- 
expansion. We obtain then the intersection numbers from this expression, 
with the scaling s 5 = ^f"S 5 , 

1 9 

< n,0 >g=l~ g, < 7"3,2 >g=2~ g 2 . 5|' 

< t- 6i o > 9=3 - <r 8i2 > g=4 - 8 5. 5! . 10? - 1^5 j 

The results up to order g = 4 had been computed in our replica article [6] 
and indeed agree with these new results, 
(iv) p > 5 

In the case p — 5 , the solutions of the conditions (12) fall in three different 
classes, (a) symmetric solution; (01,02,03,04) = (a + pi, a — pi,— a + pi,— a — 
p) with a = ±0.776887, p = ±0.321797. 

(b) cij given by ±(0.624916,-1.01437,0.53363 ±0.473515 i, 0.53363-0.473515 
i)- 

(c) dj given by ±(0.280577±0.5117 i,0.280577-0.5117 i,0.433665 ± 0.158861 
i,0.433665-0. 158861 i). 

These three cases all give a closing gap singularity with same same expo- 
nent, p(X) ~ |A| 1/5 . We have 

U(s) = J-e-T(E £> / ^ e ~ N ^ (36) 
y J Ns J 2tt v ; 

The intersection numbers may be obtained from this expression in complete 
analogy with the p=3 case. The integrand for p = 5 presents five Stokes 
lines and the spin-label j = 0, • • • , 4 characterizes the various domains with 
different asymptotic expansions. For p > 5 the situation is similar to those 
described above. 



5 Several marked points 

Up to now, we have only considered surfaces with one marked point. For 
higher intersection numbers we have to look at k-point correlation functions. 
The Fourier transform of the k-point function U(si, Sk) is also known in 
closed form ; it is given by the integral representation [14]. 

f/(si,-"-,^) = -(tre SlB ---tre SkB ) (37) 
n 

dui ^k,.„A A -A- Si s , 1 



2 Jc J /c YL 

= ( _ 1) ^-D/2 e E^ fup-jz^n na+— 3— )det 

J 1 2«7T j=lm=l M « fl m 
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For the two-point function (k=2), with the same degenerate external source 
used in (10), one has 



e N(s 1 u 1 +s 2 u 2 ) 

x- (38) 

(Ml - U 2 + Si) (Ml - U 2 ~ S 2 ) 

At the edge singularity, we expand the integrand in powers of — , with again 
the critical constraints (12), 



iV_ 

1 



p-1 

IK 1 ~ t — r) p e 



a=l °« - M 



exp[-^ T (E-)-^ T )E( ^ )] (39) 



By the shift u — > w — |s, we obtain 

^ (ai ' aa) = i72^ e 



(2m 

xe 



X (ui - «a + ^)K - «a - a±a) • (40) 

Similarly for the critical /c-point correlation function, for an arbitrary 
value of p, we have 

g _ i^l(E ^) E s » 

X "i — t w tt (41) 

det{ui — Uj + ±{Si + Sj)) 

with C = Ea=i ~ ptt • The leading connected part is obtained by the longest 
cycles (of length k) in the expansion of the determinant. 

We focus now on the case p=3 (for p—2, this integral has been already 
evaluated and the intersection numbers for two marked points for arbitrary 
genus thereby obtained in [5], (see also [10, 21]). 

For the case p=3, we consider the critical values (ai,a 2 ) = (1,-1) for 
simplicity. The singularity is located at the origin (the gap closing point). 
We replace the denominator by the integral, 

[°° dxe^ Ul ~ U2+ ^ Sl+S2)]x = \- (42) 

Jo Ui - U 2 + ^(Sl + S 2 ) 
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Then we have 



U(s u s 2 ) = / dx(e^ sl+S2 ^ x - e-^ sl+S2 ^ x ) 

Si + s 2 Jo 

f°° dVl iNsiV 3 +i{x _N s 3 )vi f°° dV2 iNs2V 3_ i{x+ N s 3 2) 



x r ClVl jNs^+Hx-^-sDv! f°° dV2 ciNs2V 3 n _ i{x+ N s 3 )v2 ^ 
J-oo 27T J-oo 27T 



which leads to 



1 f°° 1 

U{Sl > S2) = ( gl +^(^^^1/3 7o ^inh(-( 5l + 5 2 ),) 



(5i + 5 2 )(35 1 )V3(35 2 )V3 7 

where we have used the scaled variables Sj = A^ _1 / 4 §j and x = N^y. From 
now on we shall drop the tilde, but all s's should read s instead. 

Note that this expression is, from its definition, symmetric under the 
exchange of s± and s 2 although at this stage it looks asymmetrical. 

If we scale y as y — > (Ssi) 1 / 3 ^, we obtain a convenient form for expanding 
in powers of Sj, 

| roo \ 

u{s ^ ] = (., + «,)(3^,)v. y, ^(jfe+^p.o 1 ^) 

*<» - 4tW? /3 )^-0 1/3 *< - 4^"> (45) 

Note that Sj is conjugate to Aj in a Fourier transform. However, we have 
defined the dual external source model with a term Tr(BA). The definition 
(1) of the higher Kontsevich-Airy models involves a different power of A, a 
power 3 for p = 3. Thus the scaling is here given by 

Expanding then for large Aj, i.e. small Sj, the leading term is obtained by 
replacing the sinhX by its first coefficient X and dropping the s^ 3 in the 
Airy functions. Note that the corrections involve s 8 / 3 and thus are of relative 
order 1/A 8 . 

The leading term, which is order 1/A, is 

u( Sl , S2 ) = dyyAi(y)M-(-)v) 



2 s 2 Jo s 2 ' 



1 (sis 2 )»(si 3 + s 



2 si + s 2 ) 



-Ai'(0)Ai(0)) (47) 
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where Ai(0) = 3"^^ and Ai '( Q ) = _3 ^rT±)' and thus -Ai'(0)Ai(0) = |f . 
The calculation involves the differential equation Ai"(y) = yAi(y) followed 
by integrations by parts. 

The next order, which provides the intersection numbers for two marked 
points at genus one, involves three terms. The first one is the cubic correction 
X 3 / 6 to the linear term of sinh X ; the other two involve the Taylor expansion 
of the s 8//3 terms in the Airy functions. 

ac/«( Si , S2 ) = l 8l(81 + ? )2 r dyv*Mv)M-(-)*v) 

16 (3s 2 )s ■'0 s 2 



where J is 



j= r dyAi(y)M-(-)*v)- (49) 

J s 2 



AC/( 2 )( Sl , S2 ) = _^£i)i a f/ d 2 /y^'( 2/ )Ai(-(^)3 2/ ) 

O • 3 i / d S 2 S 2 

8 2 

:(^)I-%(^(0)) 2 + o L / M J (50) 



8-3 1 / sy s 2 s 1 + s 2 8-3 1 / s s 1 + s 2 

Af/( 3 )( Sl)S2 ) = -^^(-)^l(^ r dyyAi(y)Ai'(-(^)?y) 
8 • s 2 s\ Jo s 2 

- rW&S*^ - s^O* A J (51) 

The terms which are proportional to the integral J cancel beautifully 
when one adds the three correction terms. The addition leads to 

AU( Sl , s 2 ) = ^^(^2)Hsl + Sl s 2 + s 2 2 )(At'(0)) 2 (52) 

From the scaling to A, we have 

Mksl + s lS2 + 4)^^ + ^ + ^. (53) 

They correspond to £ ,o ~ X' ^ 2 >° ~ X 7 and ^.o ~ X*' and they give the 
intersection numbers for two marked points. 

< r , T 2 ,o > g =i=< Ti, >= y^- (54) 

which agrees with the solution of the Virasoro equation for p = 3 [22, 20]. 
Note that here instead of the recursive calculation of these numbers used by 
previous authors, we have, with the integral representation (45), a generating 
function for arbitrary genera of intersection numbers with two marked points. 
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6 The p ^ —1 limit 

Up to now p, which characterizes the spin-structure, was an integer larger 
than one. It is interesting to consider how the theory is modified when it is 
continued to negative values of p and in particular when p — — 1. Since the 
generalized Kontsevich model involves ^^tiB p+1 , the limit p — > — 1 gives a 
logarithmic potential trlog-B. This logarithmic potential corresponds to the 
Penner model [23, 24] and it is known to be related to the Euler character x 
of the Riemann surfaces. 

As we have discussed, the intersection numbers do depend upon p ; for 
instance, in the case of one marked point, we have [6] 

< 7i,o > fl =i= (55) 

The analytic continuation over p to negative values is thus possible and in 
this case, it gives simply < T\ q > g=1 — — ^ in agreement with the result of 

[23] : X (r\) = -T2- 

One can indeed consider continuing to negative values many of those 

formulae. For instance in our previous work concerning the replica limit, 

n — > for n x n matrix B, we had found 

limit^I < tr* >= (56, 

which is indeed finite, equal to —1/3, in the limit I — > — 1. 

The case of p — —1 is particulary interesting since Witten [4] had pointed 
out that in the limit k — > — 3 in the level-A; gauged WZW model, the inter- 
section numbers, defined by (5), become the integral of the top Chern class 
alone, 



Z 9 = -1M or 

JMg 

= -x(M Btj ) (57) 

where x is the Euler character of the manifold. 

If we compare our normalizations to Witten's correspondence with the 
level-k gauged WZW model, our definition of p is related to k by 

p = k + 2 (58) 

and the limit k — > — 3 indeed corresponds to the p = — 1 limit. 
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7 Discussion 



In this article we have shown that the intersection numbers of Riemann 
surfaces, the moduli space of surfaces with n marked points, endowed with a 
p-spin structure, are obtained from a Gaussian random matrix theory with 
external source, at an edge point where the asymptotic density of eigenvalues 
exhibits a singular behavior p(A) ~ \p . The Fourier transforms U(si, s&) 
of the k-point correlation functions provide the intersection numbers through 
the conjugacy relation Sj ~ jp. 

i 

We have found an integral representation for U(si, s^) at the critical 
point. Witten's conjecture is that the intersection numbers may be obtained 
recursively through Virasoro equations. Our formula is instead a closed ex- 
pression for arbitary p and arbitrary genus ; it also gives the possibility to 
continue in p to negative values (for instance p — — 1). It is amusing to ob- 
serve that the Fourier transform U(si, ...,Sk) has an expansion as products 
of powers of Sj, which corresponds precisely to the expansion in t n j of the 
free energy in the generalized Kontsevich model. 

Recently, this p-ih generalized Kontsevich model has been discussed as 
an effective theory of open strings between Liouville D-branes [25, 26]. The 
duality, which we have discussed with external source, corresponds to the 
relation between closed string (gravity) and open string (gauge theory) with 
cosmological constants A. The random matrix model with external source 
gives thus a theory for the case of DO branes. 

A possible extension of this work would deal with the orthogonal-symplectic 
Gaussian matrix models with external source, which are also mutually dual 
[27, 28, 29]. The time dependent case should also be reconsidered [14] at the 
light of intersection numbers theory. This is all left to further work. 
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